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Abstract 

Recently, we evaluated the virtual cross-section for Bhabha scattering in 
pure QED, up to corrections of order a 4 (Np = 1). This calculation is valid 
for arbitrary values of the squared center of mass energy s and momentum 
transfer t; the electron and positron mass m was considered a finite, non van- 
ishing quantity. In the present work, we supplement the previous calculation 
by considering the contribution of the soft photon emission diagrams to the 
differential cross-section, up to and including terms of order a^{Np = 1). 
Adding the contribution of the real corrections to the renormalized virtual 
ones, we obtain an UV and IR finite differential cross-section; we evaluate 
this quantity numerically for a significant set of values of the squared center 
of mass energy s. 
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1 Introduction 



The relevance of the Bhabha scattering process (e + e~ — > e + e~) in the study of the 
phenomenology of particle physics can hardly be overestimated. This is due to the 
fact that Bhabha scattering is the process employed to determine the luminosity 
of all the present e + — e~ colliders, at both high (~ 100 GeV) and intermediate 
(~ 1 — 10 GeV) energies, as well as in a future linear collider. For this reason, Bhabha 
scattering has been studied in great detail within the context of the electroweak 
Standard Model (we refer the interested reader to [Tj and references therein). 

In recent years, it was pointed out that, while the one- loop radiative corrections 
to Bhabha scattering and the corresponding real emission corrections were well 
known within the full electroweak Standard Model [3 E] , the two-loop corrections 
had not been calculated even in pure QED, although a large amount of work was 
devoted to the study of the contributions enhanced by factors of \n.(s/m 2 e ) A\. The 
reason for this situation was identified in the technical problem of calculating the 
necessary two-loop box diagrams. Several groups started to work on different aspects 
of the problem. In [S], the two- loop QED virtual cross-section in the limit of zero 
electron mass was calculated, while in [0] the IR divergent structure of that result was 
carefully studied. A very interesting and ambitious project aiming at the complete 
evaluation of the Bhabha scattering cross-section in QED, without neglecting the 
electron mass m, was presented in jjj. Some of the necessary Feynman diagrams 
were calculated in jHJ EH ■ 

By employing the results of jHJ E] , it was possible to complete the calculation 
of the virtual Bhabha scattering unpolarized differential cross-section, including the 
contribution of two-loop graphs involving a closed fermion loop (conventionally in- 
dicated as corrections of order a A {N F = 1), where a is the fine structure constant) 
|llj . The cross section presented in ^T] is valid for arbitrary values of the squared 
center of mass energy s and momentum transfer t. The full dependence of the cross- 
section on the electron and positron mass m was retained. In calculating the Feyn- 
man diagrams, both UV and IR divergences were regularized with the continuous 
D-dimensional regularization scheme |12j, while the diagrams were evaluated analyt- 
ically in [HI E] by means of the Laporta algorithm and the differential equations 
technique for the evaluation of the master integrals [Hj. In [TT], the renormalization 
program was carried out in the on-shell scheme; therefore the final result, expressed 
in terms of 1- and 2-dimensional harmonic polylogarithms (HPLs, 2dHPLs) ^5] of 
maximum weight 3, is free from UV divergences. However, the cross-section of ^T] 
still includes IR divergent terms that appear as a pole in (D — 4). 

As is well known, the IR divergent term in the cross-section of jTTj cancels if one 
adds to it the contribution of events of the type e + e~ — > e + e~7, in the limit in which 
the photon in the final state carries an energy which is small with respect to the 
squared center of mass energy. Such events are commonly referred to as soft photon 
events, and their contributions to the cross-section are known as real (as opposed to 
virtual) corrections. 

The purpose of the present paper is to calculate the contribution of the the soft 
radiation diagrams (up to and including order a A (Np = 1)) to the Bhabha scattering 
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differential cross-section, in order to verify that the IR divergent terms present in 
this contribution cancel against the IR divergent terms in the cross-section in jTT] . 
as well as to ultimately obtain the IR finite differential cross-section by adding the 
virtual and real corrections. 

In calculating the contribution of the soft photon emission diagrams to the cross- 
section, we integrated over the soft photon phase space, taking into consideration 
the emission of photons with an energy smaller than a certain energy threshold u. 
This threshold is supposed to be small with respect to the beam energy E. In other 
words, we assumed that the idealized detector, to be used to measure the differential 
cross-section we calculated, can tag all the events in which one photon appears in 
the final state, provided that the photon carries an energy larger than the chosen 
threshold. We also assumed that all the events in which a photon with an energy 
smaller than the threshold appears in the final state, are for the idealized detector 
indistinguishable from the Bhabha scattering events without photon emission. This 
is the standard textbook approach to the calculation of the real corrections to a 
given cross-section; admittedly, it is not enough to consider cross-sections measured 
in realistic experiments, where other aspects must be taken into account (e. g. hard 
bremsstrahlung effects, detector geometry). Very often, the experimental set up 
is so complex that the only effective tool to obtain a realistic cross-section is the 
Monte Carlo method. Nevertheless, we thought it useful to calculate the soft photon 
emission following this approach, in order to diagrammatically show how the can- 
cellation of the IR divergences works, as well as to provide a benchmark for future, 
more realistic numerical calculations. 

Instead of presenting the lengthy analytical expression of the IR finite differential 
cross-section, we preferred to implement computer codes that evaluate the differen- 
tial cross-section at order a 4 (Np = 1) for arbitrary values of the beam energy and 
scattering angle in the center of mass, E and 9, respectively. 

We have found that the corrections at order a 4 (Np = 1) are positive and very 
small with respect to the corrections of order a 3 (which are negative), so that their 
constructive contribution is strongly suppressed in the energy range of interest at 
present and future colliders. The relative weight of the corrections of order a 3 and 
a 4 (Np = 1) does increase in magnitude with the beam energy and, at a given 
energy, with the scattering angle. In order to check our numerical consistency, we 
reproduced the results in j2] and jSj revelant for our purposes. 

Our routines for the numerical evaluation of the Bhabha scattering cross- sect ion 
up to corrections of order a 4 (Np = 1), written both in Mathematica and in 
Fortran77, can be obtained from the authors |18j . 

The paper is organized as follows: in Section |2l we discuss the calculation of the 
soft real corrections to the Bhabha scattering differential cross-section 1 , up to and 
including terms of order a A (Np = 1)- In addition, we investigate the cancellation 
of the IR divergent terms of the Bhabha scattering virtual cross-section calculated 
in ^T] against the IR divergent terms originating from the soft photon emission 
diagrams considered here. In Section El we present the numerical results obtained 

^^Hcrc and in the following, we always consider differential cross-sections summed over the spins 
of the final state particles and averaged over the spin of the incoming electron and positron. 
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by evaluating the IR finite Bhabha scattering cross-section (given by the sum of 
virtual and soft corrections up to terms of order a 4 (iVp = 1)), for different significant 
choices of the beam energy E. Also, we compare the complete cross-section with 
its expansion in powers of the electron mass. We find that the first term in the 
expansion fits to sufficient approximation the numerical value of the complete cross- 
section for all the beam energies relevant in present and future colliders. In Section^ 
we present our conclusions. Finally, two appendices include the expression of the 
integrals occurring in the evaluation of the soft radiative corrections relevant to our 
calculation, as well as the leading terms in the expansion in the limit m — * of the 
Bhabha scattering differential cross-section we studied. 



2 The Real Corrections 

In this Section, we discuss the calculation of the real corrections due to the emission 
of a soft photon to the Bhabha scattering unpolarized differential cross-section in 
pure QED. In particular, we obtain the real corrections of order a 3 and a 4 (Np = 1). 
In both cases, we must consider events involving a single soft photon in the final 
state: 

e" (pi) + e + (pa) — > e" (p 3 ) + e+ (p 4 ) + 7 (*) , (1) 

where pi, p 2 , P3, p 4 , and k are the momenta carried by the incoming electron, 
incoming positron, outgoing electron, outgoing positron, and outgoing soft photon, 
respectively. All of the particles in the initial and final states are on-shell, so that 
pi = —m 2 (i = 1,4) and k 2 = 0. Therefore, we introduce the quantities s and s' 
defined as follows: 

s = -(pi+p 2 ) 2 , s' = -{pz+ptf. (2) 

By definition, the soft photon approximation consists of neglecting k in the numer- 
ator of the scattering amplitude and of setting s' = s everywhere. In this approx- 
imation, the kinematical relations that link the Mandelstam invariants s, t, and u 
to the beam energy (E) and to the scattering angle in the center of mass frame (8) , 
are 

s = -P 2 = -(pi+p 2 ) 2 ^-(p 3 + p 4 ) 2 = 4£ 2 , (3) 
t = -Q 2 = -( Pl - p 3 ) 2 ^ -(p 2 -p 4 ) 2 = -A (E 2 -m 2 ) sin 2 ^, (4) 

-V 2 = -(p! - p 4 ) 2 w -(pa - p 3 ) 2 = -4 {E 2 - m 2 ) cos 2 ~ , (5) 



u 



with 



s + t + u = Am 2 . (6) 



In the following, we often employ the dimensionless variables x, y, and z, related to 
the Mandelstam invariants s, t, and u, by the relations 



2 (l + x) 2 y/s — Vs — 4m 2 



m , x = 1^=^= , (7) 

x Vs + Vs - 4m 2 



3 



u 



-m 



— m 



:i-y) 2 

y 

:i-^) 2 

z 



\/Am 2 - t - y^t 



A/4m 2 - t + 
V4m 2 — u — y/—u 
a/ 4m 2 — u + \J — u 



(8) 

(9) 
(10) 



which are valid in the physical region s > Am 2 , t,u<0. 

The complete Bhabha scattering differential cross-section up to order a A (Np = 1) 
can be written as follows: 



da T (s,t,m 2 ) dao(s,t,m 2 ) fa\ da\ (s, t, m 2 



+ (- 

-7T 



fa\ 2 do%(s, t, m 2 



dfl dn \7T/ dfl 

where ao(s,t,m 2 ) is the tree-level (Born) cross-section 



dn 



da (s, t, m 2 
dTl 



and aT(s, t, m 2 ) 



a 



s i s 



2\ 2 



st+ — +(t- 2m 2 ) 



+ 



f 2 



2^ 2 



st + — + (s- 2m 2 ) 



+^ [(^ + t) 2 -4m 4 ]} 



(12) 



1,2) are the sum of the virtual and real corrections at order 



a and a (Np = 1), respectively. 



2.1 Real Corrections at Order ck 3 

We first consider the order a 3 contribution to the cross-section, which is given by: 
'a\ daf(s, t, m 2 ) 
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dn 



daY(s,t,m 2 ) daf(s 1 t 1 m 2 ) 



dn 



dn 



(13) 



where the superscripts V and S stand for "virtual" and "soft" . 

The one-loop virtual cross-section da\ /dVL can be found in Eq. (67) of [llj; we 
devote the remaining part of this subsection to the calculation of daf /dn. 

The diagrams contributing to the real corrections to the Bhabha scattering cross- 
section at order a 3 are shown in Fig. Q the real photon can be emitted by any of 
the incoming or outgoing fermion lines of the s- and t-channel Bhabha scattering 
tree-level diagrams. 

At this stage, it is convenient to introduce the quantity 



dermis, t, m 



dn 



da (s,t,m 2 )_^^ D ^^da^ D 4) (s,t,m 2 ) 



dn 



dn 



+ o((d-a)' 



(14) 



where: 



d(jQ (s,t,m 2 ) 

dn 



a 



r s 2i 



+ 



t 2 



1 

+ 7t 



-{s + tf - -st-m 2 {s + t) 



(15) 
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do® j dQ is the Born Bhabha scattering cross-section obtained by calculating the 
traces over the Dirac indices in D dimensions. The contribution of the s- and t- 
channel diagrams and of their interference to the r.h.s. of Eqs. ()14I15|) is explicit. It 
is then straightforward to show that, in the soft photon approximation, the contri- 
bution of the diagrams in Fig. ^ to the unpolarized differential cross-section is given 
by 

r a\ dcrf(s,t,m 2 ) /a\ da^s, t, m 2 ) 



- 



.TiJ dQ \ttJ dfl 

where the IR divergent quantity is defined as 

Jij = e { ej (pi-pj) Iij, (17) 

with ti = +1 for % — 1,4 and = —1 for i = 2,3, and where 1^ indicates the 
integral: 

1 m D - A r d D k 1 , , 



h 



lJ T (3 - f ) tt(^ 4 )/ 2 4tt 2 J k (pi ■ k) (pj ■ k) 

In Eq. f!18|) . D is the dimensional regulator; furthermore the superscript on the 
integral sign indicates that the integration should be taken over the region \k\ = 
ko < uj, with id representing the cut-off on the unobserved soft-photon energy. The 
integral in Eq. (fTSj) can be evaluated according to the standard technique discussed 
in detail 2 in Ref. |16) . It is important to observe that the integrals Iy depend only 
on the scalar product Pi ■ pj (aside from an obvious dependence on E and m), so 
that 

Hj = Iji i hi — 1 22 — I33 = 1 44 , 
Il2 = -^34 , Il3 — -^24 , Il4 = -^23 • (19) 

Consequently, the quantities Jy also satisfy the same symmetry relations. Therefore, 
Eq. (fTp'j) becomes 

(a\ daf(s,t,m 2 ) _ fa\ da°(s,t,m 2 ) ^ 

Kn) dh " U dh 4 ^ Jlj ■ (20) 

The explicit expressions of iy (j = 1,4) can be found in Appendix 1X1 

2 Thc integral in Eq. (|18|l and the one in Ref. ^Sj differ by a normalization factor. In partic- 
ular, the integral with the normalization employed in this paper can be obtained by multiplying 
Eq. (1.230) in Ref. 16 by the factor 

D-A 



where 7 is the Euler constant and \i the 't Hooft scale. Our choice of the normalization constants 
removes, in all the divergent integrals, the finite terms associated with the use of dimensional 
regularization and sets the 't Hooft scale to m; this choice is consistent with the normalization 
employed in Ref. 
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(a) 



(b) 



(c) 



(d) 




Figure 4: Diagrams contributing to the real corrections at order a 4 (Np = 1). 

The IR pole that originates from the integral I\j is multiplied, in Eq. (|20p. by 
the terms proportional to (D — 4) in o~q ; this product provides a finite contribution 
to <rf . Terms proportional to (D — 4) in erf must be neglected. 

It is useful to understand how the cancellation of the IR poles works from a dia- 
grammatic point of view. Fig. Eland Fig. El schematically describe the situation. The 
contribution to the differential cross-section of the interference of the two diagrams 
shown in the first term of each line is IR divergent. This IR divergence cancels 
against the contribution to the real cross-section given by the second term in each 
line, where the product of the two tree-level diagrams represents the contribution 
of their interference to the cross-section in Born's approximation. We remind the 
reader that the one-loop photon self-energy diagrams are IR finite. 



2.2 Real Corrections at Order a 4 (N F = 1) 

The order a 4 (Np = 1) contribution to the cross-section is given by: 



/a\ 2 da%(s, t, m 2 ) ( a \ 2 
W dTl = W 



da\ (s,t,m 2 ) dcrf (s, t, m 2 ) 



dtt 



dtt 



(21) 



where the two-loop virtual cross-section da\ /dtt can be found in Eq. (68) of 

The IR divergencies present in the Bhabha scattering differential cross- sect ion 
at order a A (Np = 1) cancel against the contribution to the real cross-section of 
the interference of the diagrams in Fig. 0] with the single photon emission tree- level 
diagrams in Fig. ^ In discussing the soft corrections to the cross-section at order 
a 4 (N F = 1) it is convenient to introduce the quantity 



derf (s, t, m 2 



dtt 



daY(s, t. 



m 



(U,S) 



dtt 



(u,s) 



+ {D-4) 



da\ 



CD-4), 



s, t, m 2 



dtt 



(U,S) 



(22) 



The first term in the r.h.s. of the equation above is the contribution to the virtual 
cross-section of the one-loop self-energy diagrams and corresponding counter-term 
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diagrams: 



da{(s, t, m 



dtt 



(u,s) 



a 



1 



st + — + (t- 2m 2 ) 2 



2R e n 



(11,0), 



t , 

st H h Is — 2m 

2 v 



2\2 



2n 



(11,0) 



(t) 



+1 [(* + t) 2 - 4m 4 ] (Reltf'' ^) + D^t)) }, (23) 



where IIq 1/ ' ' ) is the UV renormalized photon self-energy. The latter quantity has 
been discussed in detail in Section 4 of ^T]. The second term in the r.h.s. of 
Eq. ((22) is: 



da[ D 4 \s,t,m 2 



dQ 



(u,S) 



a 



1 



st + — + (t- 2m 2 ) 2 



2Rer4 1M) ( S ) 



t 2 



st H h (s — 2m 

2 y 



2\2 



2R 



(*) 

J-[(s + t) 2 -4m 4 ] (Ren ^( S ) + nr)(t) 



si 

Ren ( 



8 ) + n o"'° ) ( t ) +![(, + t) 2 - s t-2m 2 (,+t)] x 



x (Ren^ 0) ( s ) + n^ 0) (t) 



(24) 



where ELq j1 is the term proportional to (D — 4) in the expantion of the renormalized 
photon self-energy (the explicit expression of this quantity can be found in the 
appendix of [TT]). 

The contribution of the soft photon emission to the real corrections to the Bhabha 
scattering cross-section at order a A (N F = 1) is given by 



a\ 2 dcr|(s, t, m 2 ) ( a\ 2 ( do^ '(s,t,m 2 ) 

TV 



dn 



7T 



dn 



(li,S) 



3=1 



lj > 



(25) 



where the integrals J\j (j — 1, • • • , 4) have been introduced in the previous subsec- 
tion. The term proportional to (D — 4) in Eq. (|2*2*|) provides a finite contribution to 
the real corrections in Eq. (|25p. since Jy contains an IR pole. Terms proportional 
to (.D — 4) in erf (s, t,m 2 ) are then neglected. 

Figs. |3HH1 explicitly show how the cancellation of the IR divergences takes place 
from a diagrammatic point of view: the contribution to the virtual cross-section of 
the interference of the diagrams in the first term of each line is IR divergent; such 
divergence cancels against the interference of the two diagrams in the second term 
multiplied by the appropriate combination of J\j integrals. We observe that in the 
last two lines of Fig. in the last line of Fig. and in the second line of Fig. [3 
the subtraction of the real radiation in the second term of the l.h.s. does not cancel 
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Figure 5: Cancellation of the IR divergencies of the two-loop Np = 1 reducible 
diagrams. 
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Figure 6: Cancellation of the IR divergencies of the products of one-loop self-energy 
and vertex diagrams. 
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Figure 9: Differential cross-section for E = 22 GeV. The dashed- dotted line corre- 
sponds to the Born cross-section, the dashed one to the cross-section up to correc- 
tions of order a 3 , and the continuous one to the cross-section up to and including 
corrections of order a 4 (Np = 1). 



completely the IR pole in the corresponding virtual correction (first term in the 
l.h.s.). A residual IR pole, proportional to C(2), remains. As expected, the sum of 
the residual poles vanishes and the cross-section is therefore IR finite. 

3 Numerical Results 

In order to numerically evaluate the Bhabha scattering cross-section up to correc- 
tions of order a 4 (N F = 1), we developed two computer codes. One of them was 
written in Mathematica |T7j, while the other was written in Fortran77. Follow- 
ing [3], in the numerical calculation, we fixed the energy cut-off on the undetected 
soft photon to u = 0.1 E, where E is the beam energy. We compared the results of 
the two codes for a beam energy ranging from E = 0.01 GeV to E — 500 GeV, and 
for an arbitrary choice of the scattering angle 9, finding complete agreement. 

In Fig. |UJ we show the differential cross-section for E = 22 GeV and for < 6 < 
7i. The dashed-dotted line represents the cross-section in the Born approximation 
(Eq. (|12)l). while the continuous (dashed) line represents the cross-section at order 
a 3 (a 4 (ATp = 1)), Eqs. (jl3!21|) . The radiative corrections at order a 3 are negative; 
therefore, they lower the cross-section. The corrections at order a A (Np = 1) are 
negative and very small with respect to the corrections of order a 3 ; they lower 
the differential cross-section, although, as seen in Fig. the effect is difficult to 
appreciate graphically. 

The relative weight of the corrections of order a 3 and a 4 (Np = 1) are shown 
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Figure 10: Ratio Ri of the one-loop corrections over the tree-level cross-section. 
The different lines correspond to the following beam energies (from top to bottom): 
E = 10 MeV, 100 MeV ; 1 GeV, 22 GeV ; 100 GeV ; and 500 GeV. 

in Fig. El and Fig. ^2 for six different choices of the beam energy. In Fig. El we 
plotted the ratio of the order a 3 corrections to the tree-level cross-section: 



It is evident that the relative weight of the correction increases in magnitude with 
the beam energy and, at a given energy, with the scattering angle. A similar plot 
can be found in |3J, for E = 22 GeV, where the full set of one-loop corrections in 
the Standard Model were considered. At 9 = n/2 the corrections range from —9% 
of the Born cross-section at E = 10 MeV to -37% for E = 500 GeV. 

In Fig. HU we plotted the ratio of the order a A (N F = 1) corrections to the 
complete order a 3 cross-section: 



The corrections increase in magnitude with the energy and, for fixed energy, with 
the scattering angle. The relative weight of these corrections at 9 = ir/2 ranges 
from —0.08% of the complete cross- section at order a 3 for E = 10 MeV to —4.3% 
for E = 500 GeV. 

Figs. ^1 and El show the dependence of the cross-section on the beam energy, 
for small and large scattering angles respectively. Similar plots for a wider choices 
of angles can be found, limited to the one-loop corrections, in [2]. Using our codes, 
we reproduced the plots shown in j2] finding agreement. 




(26) 




(27) 
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Figure 11: Ratio R2 of the two-loop corrections over the tree-level cross-section. The 
different lines correspond to the following beam energies (from top to bottom): E = 
10MeV ; 100 MeV ; 1 GeV, 22 GeV, 100 GeV ; and 500 GeV. 




50 100 150 200 250 300 350 400 450 500 

E (GeV) 

Figure 12: Bhabha scattering cross-section as a function of the beam energy for 
a scattering angle 6 = 5°. The dashed-dotted line corresponds to the Born cross- 
section, the dashed one to the cross-section up to corrections of order a 3 , and the 
continuous one to the cross-section up to corrections of order a 4 (Np = 1). 
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Figure 13: Bhabha scattering cross-section as a function of the beam energy for a 
scattering angle 9 = 90°. The dashed- dotted line corresponds to the Born cross- 
section, the dashed one to the cross-section up to corrections of order a 3 , and the 
continuous one to the cross-section up to corrections of order a 4 (Np = 1). 



To complete the analysis of the numerical results, we expanded the analytic 
expression of the cross-section in the limit in which the squared electron mass is 
negligible with respect to the kinematic invariants s, t, and u. We define the leading 
terms of the cross-section in the limit m 2 — > through the relations: 



daj 



T 



da] 



dtt dtt 



+ o 



m m m 
s t u 



f28) 



where % = 1,2 and where the subscript "L" stands for "leading". The expressions 
of (daf /dfl)\L can be found in Appendix iBl 

In Figs. HU and [T3] we plotted (for a fixed value of the scattering angle) the 
quantities: 



Di 
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dQ 
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dn 



dn 



da, 
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o a da 1 
dQ 7T dQ 



-i 



(29) 
(30) 



respectively. A glance of the figures shows that the leading terms of the cross-section 
approximate to a continuosly better degree the complete results for increasing beam 
energy. The leading terms of the cross-section fail to reproduce the complete result 
in the extremely forward and backward regions, where t or u becomes smaller than 



m 
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Figure 14: D\ for 9 = 5° (dashed line) and 90° (continuous line). 
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Figure 15: D2 for 8 = 5° (dashed line) and 90° (continuous line). 
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The codes for the numerical evaluation of the Bhabha scattering differential 
cross-section in pure QED up to order a 4 (Np = 1) are available from the authors 

m 

4 Conclusions 

In the present paper, we completed the evaluation of the Bhabha scattering cross- 
section at order a 4 (Np = 1) in pure QED. The calculation was performed without 
neglecting the electron mass m, and is valid for all physical values of the independent 
Mandelstam invariants s and t. The master integrals necessary for the evaluation 
of the virtual corrections were calculated in [Hj and [Oj, while the UV renormal- 
ized unpolarized differential cross-section was obtained in jTIJ. The calculation was 
completed by providing the real correction in the approximation of a soft photon 
emission up to order a A (Np = 1), as well as by explicitly showing that the IR poles 
present in such corrections cancel the remaining IR poles of the virtual cross-section 
calculated in [TI]. Finally, we developed computer codes for the numerical eval- 
uation of the UV and IR finite cross-section. We compared our findings with the 
results present in the literature (where possible), obtaining complete agreement. We 
verified that the effect of the terms proportional to positive powers of the electron 
mass m is negligible in the energy range of interest in present and future colliders. 
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A Integral Iij 

As explained in Section EJ the integrals iy are symmetric with respect to their 
indices and satisfy the symmetry properties listed in Eqs. (|19|) . Therefore, in order 
to calculate the real corrections to the Bhabha scattering cross-section up to the 
order that we consider in this work, one needs to calculate the integrals I\j (j = 1, 4). 
Following the procedure outlined in Ref. ^B] , it is straightforward to verify that, 
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after the integration on the photon phase space has been carried out, the integral 
becomes 



Pij 



1 



1. (u 2 



D-4 2 



In 



<->) /[f-A/,, 



+ £>(£> -4). (31) 



In the equation above, the quantity pjj has a particularly simple expression in terms 
of the dimensionless quantities x, y, and z introduced in Eqs. (JZHHl): 



111 

Pl2 = - , Pl3 = - , Pl4 = - , Pll = 1 • 

X y Z 



(32) 



The quantity in Eq. (|31|) can be expressed as a simple integral as follows: 
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After observing that the scalar products p\ ■ pj also have a simple form in terms of 
x, y, and z, 
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In addition, the quantity A/y can also be expressed in integral form: 



A J 
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where the Lorentz vector P^- is defined by the relation 
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The integral in Eq. (J3~9")) can also be evaluated according to the procedure outlined 
in Ref. ^B]; one finds: 
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where I = 2, 3, 4 and the arguments af^ (k = 1, . . . , 4) have convenient expressions 
in terms of x and pu : 
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,( 4 ) 
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We observe that and are equal to zero; as a consequence, the corresponding 
dilogarithms in Eq. (j42j) also vanish. 



B Leading Corrections at Order a 3 and Order 
a 4 (N F = 1) 

In this Appendix, we provide the explicit expressions of the leading radiative cor- 
rections defined in Eq. (J2HJ) - 
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